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Abstract 

We numerically investigate the decay of initial quantum Fock states and their superpositions for 
a mechanical resonator mode coupled to an environment comprising interacting, damped tunneling 
two level system (TLS) defects. The cases of one, three, and six near resonant, interacting TLS's 
are considered in turn and it is found that the resonator displays Ohmic bath like decay behavior 
with as few as three TLS's. 

PACS numbers: 85.85.+j,03.65.Yz 



I. INTRODUCTION 



The quest to understand the quantum-to-classical transition has led to the development 
of macroscopic mechanical systems in which researchers hope to realize quantum states. 
In a 2010 landmark experiment,^ a state corresponding to a single quantum of vibrational 
energy in a mechanical resonator was created and its subsequent decay dynamics measured. 
We anticipate that similar measurements involving higher number quantum Fock states in a 
mechanical system will be achieved in the near future. In light of these developments, there 
is a need to understand the decoherence mechanisms in these systems. 

In the early 70 's, researchers studying the unusual thermal and acoustic properties of 
amorphous insulators at cryogenic temperatures suggested that the materials might be pop- 
ulated by tunneling two-level system (TLS) defects.^ More recently, evidence for the rele- 
vance of TLS's in micronscale superconducting qubit dynamics has been established.™ The 
same amorphous materials are often used in the fabrication of micronscale mechanical sys- 
tems and thus it is likely that TLS's will play a dominant role in their quantum-to-classical 
transition at cryogenic temperatures. 

In particular, we anticipate that TLS's will provide one of the main mechanisms for the 
decay of quantum states in mechanical resonators. In Ref. [12J, we presented an estimate 
indicating that a given low order flexural mode of a micronscale mechanical resonator vi- 
brating at radio frequencies may be near resonance with a few TLS's, but is unlikely to 
interact resonantly with large numbers of TLS's. These TLS's couple to the motion of 
the resonator via its strain, and thus will be part of the environment responsible for the 
decay of quantum flexural modes. Reference [12] numerically investigated the damping of 
initially coherent states and the decoherence dynamics of initial superpositions of spatially 
separated coherent states, where the environment consisted of either one or three damped 
TLS's. Clear signatures of resonator amplitude dependence were observed in the damping 
and decoherence dynamics, a consequence of TLS saturation. This behavior is qualitatively 
different from the amplitude-independent damping and decoherence resulting from the stan- 
dard, Ohmic oscillator bath model of an environment. However, it is of interest to explore 
the damping and decoherence dynamics as the number of near resonant TLS's increases, in 
particular to establish the expected transition to approximately Ohmic like behavior. 

In this work, we numerically model the low temperature (hu> <C ksT) damping and 
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decoherence dynamics of a mechanical resonator coupled to between one and six damped 
TLS's that are near-resonant with the resonator, where the latter is initially prepared in 
either a single Fock state or superposition of Fock states. We find, perhaps surprisingly, that 
the damping and decoherence dynamics resembles quite closely that for the Ohmic oscillator 
bath model even with only three near resonant damped TLS's furnishing the mechanical 
resonator environment. In particular, the Fock state lifetime is observed to scale closely 
as 1/n, where n is the initial number of resonator quanta (Fock state number), while the 
decoherence time of a superposition of ground and excited Fock states is found to be close to 
twice the decay time of the excited state, both in accord with the Ohmic model. A partial 
understanding of these numerical results can be obtained from a simpler, Born-Markov 
approximated master equation model for the resonator subsystem that treats perturbatively 
the coupling between the resonator and damped TLS's to second order (with the latter traced 
over as the bath) and which facilitates analytical calculations for the decay times. However, 
even more surprising is the observation that completely removing the TLS's damping does 
not alter the Ohmic trends, even though the Born-Markov master equation model is no 
longer valid. The latter observation is reminiscent of recent numerical investigations to 
establish subsystem thermalization of closed, many-body quantum systems.^ 

In the next section we present our model system-environment master equation, with a 



more detailed derivation given in Ref. [12]. Sec. Ill investigates the damping dynamics of 
Fock states and decoherence dynamics of superpositions of Fock states for a mechanical 
resonator coupled to first a single TLS, then three near-resonant TLS's, and finally six near- 
resonant TLS's, where direct interactions between the TLS's are neglected. An approximate 
master equation model is presented, yielding analytical decay rate expressions that partially 



explain the numerically observed trends. In Sec. |IV[ we begin with deriving the oscillator- 
TLS Hamiltonian with pairwise interactions between TLS's taken into account. The effect 
of TLS-TLS interactions on the damping of Fock states and decoherence of Fock state 
superpositions in a resonator coupled to first three and then six near-resonant, interacting 
TLS's is investigated. Finally, we offer some concluding remarks in Sec. |Vj 
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II. RESONATOR-TLS SYSTEM EQUATIONS 



In this section we present the model for the resonator-TLS system. For the TLS Hamil- 
tonian we have 



#tls — ^ 
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(1) 



where a = 1,2, ...,N labels the TLS, Aq is the asymmetry of the ath TLS's potential well 
and Aj is its tunnel splitting that depends on the well barrier height and width. Writing 
out the full oscillator-TLS system Hamiltonian, we have 

N 



H s = hu{a^a + 1/2) + 



a=l 



(2) 



where a) and a are raising and lowering operators for the resonator mode of interest, sat- 
isfying the commutation relation [a, a)] = 1. The strength of the coupling depends on 
the location of the TLS defect within the resonator. 

In Ref. [12] we derive the following master equation describing the dissipative dynamics 
of the coupled resonator-TLS system: 

Ps® = -i^^-S [y ' {Pr ' Ps(0}] "^ C ° th (2^) 
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where ps{t) is the resonator-TLS system density matrix, Y = Y zp (a+a*) gives the mechanical 
resonator mode displacement, with Y zp the zero-point displacement uncertainty, Py is the 
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resonator mode momentum, and E^ = y (Aq°^) 2 + (A^) 2 is the ath TLS energy level 
separation. The parameter 7 gives the energy damping rate of the oscillator in the absence 
of the TLS's, while the parameter T± a ^ gives the ath TLS relaxation time from its excited 
energy eigenstate in the absence of the oscillator. We shall use dimensionless time units, 
t — > out, with T\ and 7 expressed as uT\ and j/u, respectively, and A, A», and temperature 
T expressed as Xjfko, Ai/ftw, and hsT/hco, respectively. 
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III. DAMPING AND DECOHERENCE DUE TO NON-INTERACTING TLS'S 



A. Single TLS 

In this section we investigate the damping of Fock states and the decoherence of Fock 
state superpositions in a mechanical resonator interacting with a single damped TLS. As a 
partial check of our numerical methods, we begin by evaluating the number state probability 
Pn — ( n \p\ n ) as a function of time for a resonator mode coupled to an Ohmic oscillator bath 
only, where the analytical solution is known. Fig. [T] shows the log of the number state 
probability for initial Fock states \n), with n = 0ton = ll, in the absence of the damped 
TLS. The slope of each successive curve decreases by an increment of 1/Tn, where Tn is the 
lifetime of the first excited state; as expected, the number state lifetime decays as 1/n.™ 
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FIG. 1: (Color online) Log of number state probability P n vs tot for a range initial Fock states for 
the resonator coupled to an Ohmic oscillator bath only, where 7 = 0.01 and T = 0.09. 

Fig. [2] shows the number state probability for the resonator coupled to a single, on- 
resonance damped TLS. In this case we see that the P n curves oscillate, as energy is trans- 
ferred from the resonator to the TLS and back. As a partial check of the numerics, the 
time of the first minimum of each curve corresponds closely to the Jaynes-Cummings model 
prediction for the transfer time of a quantum of vibrational energy to a symmetric, on res- 
onance TLS: ut = nE/(2Xy/n), where hu = E = A&. The left-hand plot shows the number 
state probability for four low-n states, while the right-hand plot displays curves for higher 
energy states. The three high-n states all appear to decay at the same rate, as indicated 
by the black curve, which simply follows the maxima of the undulating curves. The same 
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FIG. 2: (Color online) Number state probability P n vs ut for various initial Fock states for the 
resonator coupled to a damped TLS only. The black curve, indicating the peaks of the curves in 
the right-hand plot, is the same in both plots. For both plots Ao = 0, Af, = 1, A = 0.1, T = 0.09, 
and Ti = 10. 

black curve is shown in the left-hand plot, and in this case the P n plots clearly fall short 
of this "maximum" curve. Further, the extent to which the curves fall short depends on 
their n value, with the n = 1 curve having the lowest amplitude. Thus, for low-energy 
Fock states the Fock state probability appears to exhibit n-dependent damping, while for 
higher-n states the damping does not depend on the number state. This n-dependence may 
be an indication of TLS saturation. For low-ra Fock states the unsaturated TLS pulls energy 
from the resonator, where it is then dissipated to the TLS bath, causing the number state 
probability to decay more quickly. For higher-ra states, however, the TLS is saturated and 
thus contributes uniformly to damping, independent of the initial n. 

Next, we investigate the effect of a damped TLS on number state superpositions. A useful 
'visual' representation of the state is its Wigner function; Fig. [3] shows two initial oscillator 
states: an equal mixture of the ground and n = 7 state [Fig. |3(a)| , and a superposition of 
the same two Fock states [Fig. 3(b)] . In both cases the Wigner function has positive and 



negative values, because both the Fock state mixture and the superposition are non-classical 
states. However, the spoke-like interference fringes in the superposition plot indicate the 
presence of non-zero off-diagonal terms of the density matrix, as opposed to the concentric 
undulations in the mixture plot. Fig. [4] shows four equally-spaced snapshots of the Wigner 



function for a resonator initially in the superposition state shown in Fig. 3(b) The resonator 
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(a) (b) 

FIG. 3: (Color online) Wigner function for a mixture of number states |0) and |7) (a) and for a 
superposition of the same states (b) . 



is coupled to an Ohmic bath that causes the state to decay and the amplitude of the Wigner 
function to decrease. However, both the ring and spoke-like structures of the initial state 
are still visible in the final snapshot. Fig. [5] shows a similar set of snapshots, this time for a 
resonator coupled to an on-resonance, damped TLS only. In contrast to the superposition 
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FIG. 4: (Color online) Evolving Wigner function for the resonator initially in the superposition 
state shown in Fig. 3 (b)| coupled to an Ohmic oscillator bath with 7 = 0.01 and T = 0.09. 



state decay in the Ohmic bath case, we see that the spoke-like structure disappears first, 



leaving concentric rings similar to those seen in Fig. 3(a) The dephasing time T$ is usually 



defined in terms of the decay times of the on- and off-diagonal terms of the resonator's 
density matrix as follows: 

1 1 1 
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where T 0n is the lifetime of the off-diagonal density matrix element po n , and T nn is the 
lifetime of the diagonal matrix element p nn . The disappearance of the spokes prior to the 





FIG. 5: (Color online) Evolving Wigner function for the resonator initially in the superposition state 
shown in Fig. |3(b) coupled to a damped TLS with A = 0.1, Aq = 0, A& = 1,T = 0.09, andTi = 10. 



rings suggests a finite T^, in contrast to an oscillator bath, where T 0n = 2T nn . 



B. Three TLS's 

We now increase the number of damped TLS's to three. The TLS energies A ( a) and 
A[ Q \ a = 1,2,3, are chosen randomly according to the Standard Tunneling Model (STM) 
distribution] 12 * 15 ! As our condition for near resonance, the corresponding TLS energies 
are restricted to the range 0.75hu < £ {q) < 1.25fkj, where recall = yj '(A^ ] ) 2 + (A[ a) ) 2 . 
We also choose random values for the relaxation times of each individual TLS by first 
selecting a reference T\ value and then assigning to each TLS a randomly-generated 
within ±50% of the reference value. Furthermore, each TLS is assigned a random A^ 
coupling that is within ±50% of a reference value A = 0.1/6, scaled down from the single 
TLS coupling considered in the previous section (A = 0.1) so as to avoid significant TLS- 
induced renormalizations of the resonator's harmonic potential resulting from having more 
coupled TLS's. We choose a temperature T = 0.09 for all plots. 

To investigate Fock state decay, we choose an initial Fock state \ipo) = \n) and then 
determine the corresponding number state probability P n as a function of time. Fig. [6] 
shows the decay of P n for a range of initial Fock states. In contrast to the single-TLS case, 
the number state probabilities do not show large oscillations but instead decay relatively 
smoothly. Furthermore, the nearly- linear curves in the log plot indicate that P n decays 
exponentially and the decay rates can be extracted from a linear fit. We noted in the 
previous section that for a resonator damped by an Ohmic bath, the decay time for the n th 
state goes as T nn = Tu/n, where T\\ is the decay time for the first excited state: the decay 
rate scales as n. Fig. [7] shows the normalized decay rate Tn/T nn for P n as a function of n . 
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FIG. 6: (Color online) Left: P n vs Lot for a resonator coupled to three non-interacting TLS's. Right: 
Log of P n vs Lot. For all curves T\ = 10. 




FIG. 7: (Color online) Normalized decay rate vs n for single Fock states. The resonator is coupled 
to three non-interacting TLS's (solid) and to an Ohmic bath without any TLS's present (dot-dash). 
For both curves T\ = 10. 

For a resonator coupled solely to an Ohmic bath, the curve has a slope equal to one. For a 
resonator coupled to three TLS's, the slope is very close to one; Fock states decay similarly 
to a resonator that is Ohmically coupled to a bath of free oscillators. 

We now investigate the decay of a superposition of the ground state and the nth excited 
state, \ip) = 1/V2(|0) + |n», with each TLS initially in a thermal state. We consider the 
p nn and pon elements of the density matrix as a function of time (plots not shown). The 
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curves decay approximately exponentially, and the log plots decay linearly. We thus apply a 
linear fit to the natural log of the curves to find the diagonal and off-diagonal decay times, 
T nn and T 0n , respectively. Fig. [8] shows the log of the decay times as a function of log(n) for 
a range of T\ values. We plot 2T nn to allow for a comparison to the relation T$ n = 2T nn for 
an Ohmic bath. The curves in Fig. [8] all decay uniformly and with a slope ~ — 1. The 2T nn 
and T 0n curves are very similar: dephasing is negligible compared to decay. 

The curves in Fig. ^ show a surprising dependence on T\. As a reminder, is the decay 
time of the ath TLS from its excited to ground state. Because T\ determines the strength 
of the coupling between a TLS and its bath, with smaller T\ corresponding to stronger 
coupling, we would expect T nn and T 0n to decrease as T\ decreases; stronger coupling would 
result in shorter resonator Fock state decay times. However, Fig. [8] shows that the opposite 
is true. The curve with T\ — 1 shows longer decay times than the curve with T\ = 100. The 
lowest (solid green) curve is for a resonator coupled to three undamped TLS's, and thus a 
T\ for this curve is not given. This curve shows the shortest decay times, and appears to be 
the large- T\ limit of the curves for the damped TLS's. 
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FIG. 8: (Color online) To n and 2T nn for the initial superposition state |V>o) = Vv2(|0) + \n)) 
with a range of T\ values shown in the legend. For all curves the resonator is coupled to three 
non-interacting TLS's. For the lowest curve (solid green with star markers) the TLS's are not 
damped, and thus a T\ time is not given. 

To further investigate the 7\ dependence of the decay times, in Fig. [9] we plot T 0n and 

10 



600 



500 



CD 

E 400 
>, 

CO 

300 
Q 



200 



100 



. 2*T no interactions 



10 20 30 40 50 60 70 80 90 100 



FIG. 9: 2T nn and T$ n vs T\ for the n = 4 superposition state. The resonator is coupled to three 
non-interacting TLS's. 

2T nn as a function of Ti for the n = 4 superposition state. The plot shows a strong depen- 
dence on Ti, particularly for T\ < 10, and suggests that reducing the TLS-bath coupling 
causes superposition states to decay more quickly. This surprising dependence on 7\ will be 
discussed in further detail below in Sec. HI D| Finally, in Fig. 10 we plot 2T nn and To„ vs n 
for three different realizations of the randomized TLS parameters. While the curves indicate 
the same qualitative linear dependence T nn = Tu/n, there is some scatter in the T n values, 
as indicated by the different intercepts. This is to be expected given that we have only a 
small statistical sample of three randomly selected TLS's coupled to the resonator. 



Six TLS's 



We now consider a resonator coupled to six non-interacting TLS's. We assign random 
values to the TLS energies Aq°^ and according to the STM distribution, as well as 
random values to resonator-TLS coupling term and the TLS times, selected as in 
the previous section. The temperature T = 0.09 for all plots. We first consider the decay 
of a Fock state as a function of time for a resonator coupled to six damped spins. From the 



log plot in Fig. 11, we see that the natural log of the number state probability P n decays 



approximately linearly with time. The oscillations at long times for the higher energy states 
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FIG. 10: (Color online) 2T nn and To n vs n for three different groups of the TLS parameters. For 
all curves T\ = 10. 




FIG. 11: (Color online) Left: P n vs ojt for a resonator coupled to six non-interacting TLS's. Right: 
Log of P n vs uit. For all curves T\ = 10. 

are numerical artifacts arising from the exponentially small P n values. We can apply a linear 
fit to the log plot to determine the n-dependence of the decay rate. As for the resonator 
coupled to three TLS's, we find that the resonator's normalized decay rate scales with initial 



Fock state number similarly to that of an Ohmic bath, i.e. with slope ~ —1 (see Fig. 22) 



Next, we study the decay of a superposition of the ground state and the nth excited 



state. Fig. 12 shows the log of T 0n and 2T nn vs the log of the initial n characterizing the 



superposition state, for seven different values of the average TLS T\ time. We note that 
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FIG. 12: (Color online) To n and 2T nn for the initial superposition state IV'o) = l/v2(|0) + |n)) 
with a range of T\ values. For all curves the resonator is coupled to six non-interacting TLS's. For 
the lowest curve (solid blue with star markers) the TLS's are not damped, and thus a T\ time is 
not given. 



all of the curves have a slope ~ —1. Similar to Fig. Miov three TLS's, Fig. 12 shows little 



difference between T 0n and 2T nn for the different values of 7\; dephasing is negligible. 



In Fig. 13 we show the 7\ dependence of the on- and off-diagonal decay times for the 
n = 4 superposition state. The plot shows the same strong dependence on T x as for the case 
of three TLS's (Fig. [9]). As a reminder, T\ is the time it takes for a TLS in its excited state 
to decay to its ground state. Thus, we would expect that as we decreased T\ the resonator 
states would damp more quickly, resulting in a shorter decay time. For six TLS's, however, 
we find that as we decrease T 1; the diagonal and off-diagonal terms of the density matrix 
decay more slowly. This unexpected behavior suggests that the coupling between the TLS's 
and their individual baths is somehow obstructing a more efficient means of dissipation. This 



is supported by the lowest curve in Fig. 12, which is for a resonator coupled to six TLS's that 



are not coupled to their individual baths, and yet indicates the shortest oscillator Fock state 



decay time. In Sec. |III D[ we show through an analytical approximation that this behavior 
can be partially explained by considering the TLS bath Lorentzian line width dependencies 
on T\ . 



Finally, as we did for three non-interacting TLS's, we now plot the decay of T nn and T ( 



Of i 



13 



200 




10 20 30 40 50 60 70 80 90 100 



FIG. 13: 2T nn and Tq„ vs T\ for the n = 4 superposition state. 



for three different realizations of the TLS parameters. Fig. 14 shows that the groups of six 
TLS's exhibit a higher degree of agreement than the three-TLS groups did, with uniform 
slopes ~ —1. This is a good indication that we have moved to a regime more akin to a 
dense TLS spectrum, with variations in the parameters of individual TLS's having less of 
an impact on the resonator. 
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FIG. 14: (Color online) 2T nn and T$ n vs n for three different realizations of the TLS parameters. 
For all curves T\ = 10. 
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D. Analytical Approximation to Fock State Damping 



In this section we present an analysis of Fock state damping due to TLS's. We assume 
that the coupling between the mechanical resonator and the N TLS's is sufficiently weak 
that we can make a self-consistent Born approximation, where we expand perturbatively to 
second order in the resonator-TLS couplings and trace over the TLS's to obtain the following 
resonator master equation: 

Pm(t) = -- [H m ,p m (t)} 

i r* ... n 



,, 2 j dt' {-({B(t),B(t')})[Y,[Y(t' -t),p m (t)}] 
+ l([B(t),B(t')])[Y,{Y(t' - t),p m (t)}]y (5) 

where p m and H m are the mechanical resonator density matrix and Hamiltonian, respectively, 
and 

N 

B(t) = Y,^ a) 4 a) (t), (6) 

a=l 

with \( a ' the coupling between the oscillator and the ath TLS. Solving for the TLS- 
environment dynamics in the absence of the resonator, one can find the symmetric 
({B(t), B(t')}) and antisymmetric ([B(t), B(t' )]) correlation functions of the TLS bath. 
Thus, in the above Born approximation, we neglect the influence of the resonator on the 
TLS dynamics. More specifically, the approximation does not account for possible nonlin- 
ear, resonator amplitude-dependent saturation effects, or the possibility of coherent energy 
exchange between the resonator and the TLS's. The importance of these effects depends 
on the relative coupling strengths between the mechanical resonator and the TLS's, and 
between the TLS's and their respective baths. Following the analysis in Ref. j9], we have 
for the TLS bath correlation functions: 

1 N 

U{B(t),B(t')}) = ^(AWf ^^tl-^c-^ 

a=l 

-.(«)/■ 



and 



+ sin 2 cos [E^(t - t')/h] (7) 



i N 

U[B(t),B{t'))} = (A (a) ) 2 sin 2 ^(a^)sin [E^(t-t')/h] e'^^, (8) 



a=l 
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where (a ( z a) ) = tanh (E^/(k B T)), sm6^ = A { b a) /E^\ and cos#( Q ) = A { a) /E^. The TLS 
dephasing rate is given in terms of the relaxation rate as 

2" 



.(«) 
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(a) 



r (a) 
1 1 > 



(9) 



where = T[ 



We now substitute Eqs. (7) and (8h into the mechanical resonator 



master equation ([5]), and insert the free resonator (oscillator) dynamics solution 

Y(t' -t) = Y zp (ae- iuJm(t '- t) + a ] e iUm{t ' - t] ) . (10) 



We then make a rotating wave and a Markov approximation, and assume temperatures 
ksT <C such that (c^) ~ 1. We thus obtain the probability that the mechanical 



resonator is in the nth Fock state, P n 



dP n (t) 
dt 



[n\p m \n): 
7Fock[riP n (0 - (n + l)P n+1 (t)} 



(11) 



where 7Fock(= T n ) gives the decay rate for an initial n — 1 Fock state: 



7Fock 



1 N 



sin 2 ^ Q )- 



2T 



(a) 



a=l 
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+ (E( a )/h-u„ 



(12) 



Eq. (11 ) shows that the decay rate for an initial n Fock state scales with n, as we saw in the 



numerical simulations. The dependence of the probability decay rate on T\ comes from the 
TLS dephasing rate dependence of the Lorentzian term. From Eq. (9), we see that 
scales as T^. We now consider the form of the given Lorentzian, subject to the rescaling 

2T ^ (13) 



(r 2 /ef + (E/h-u m ) 2 ' 



Fig. [15] shows the Lorentzian factor as a function of u m for three different e values. As we 
increase e (i.e., increase the TLS damping time, T\), the Lorentzian factor correspondingly 
increases, as long as \E/h — u m \ < T 2 /e, i.e., within the Lorentzian linewidth. Physically, 



Eq. ( 12 ) indicates that for a mechanical resonator that is approximately resonant with a TLS, 



the longer the TLS decay time, the more rapidly it absorbs energy from the mechanical 
resonator, and hence the shorter the Fock state probability decay time. However, as e 
continues to increase, we eventually have that \E/h — u m \ > T 2 /e. The TLS is no longer 
approximately resonant with the oscillator, and so the Lorentzian factor and thus the decay 
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FIG. 15: (Color online) Lorentzian function vs uo m for three different e values. TLS parameters are 
A = 0.6281, A b = 0.7592, and T x = 0.6396. 



rate decreases. Fig. 16 shows the dependence of the decay time on e that follows from one 
of the distributions of TLS-oscillator coupling and parameter values used in the numerical 
simulations. The intermediate dip is due to some of the TLS's going out of resonance. While 




FIG. 16: Fock state decay time vs TLS T\ scaling factor e for the n = 1 Fock state. 

the plot does not show quite the same monotonically decreasing decay time with increasing e 
as found in the numerical simulation, it does give approximately the same overall decreasing 
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trend. Differences are due to the breakdown of the Born-Markov approximation for treating 
the TLS subsystem bath. 



IV. TLS-TLS INTERACTIONS 



A. Derivation of Hamiltonian 



Experiments have shown that interactions between TLS's play an important role in dissi- 
pation and decoherenceP^H^In this section we derive the TLS-TLS interaction Hamiltonian. 
We begin with the Hamiltonian for an elastic wave system interacting with TLS defects. The 
Lagrangian for an elastic wave system is^ 

L wave = i J d 3 r[pui(f,t)ui(f,t) - Cij k idiUj{f,t)d k ui{f,t)\, (14) 

where V is the system volume, p is the mass density, Ui(r,t),i = 1, 2, 3 is the ith component 
of the displacement vector field, and is the elastic modulus tensor. We use the Einstein 
summation convention. The Hamiltonian is by definition 

tt • "^-"wave j (ir\ 

-"wave -^wavej V / 

dm 



which, with Eq. [T4j gives 

(16) 



-"wave ' ' 



1 1 ' ' ' 

V 



^UiUi + ^CijkidiUjdkUi 



In addition to the non-interacting TLS Hamiltonian ([I]), we have the TLS-wave system 
interaction Hamiltonian 

N 

H ** = - E K-W a) ] > ( 17 ) 



a=l 

where v\j is the deformation potential tensor at the a TLS location and 

4 Q) = \[9iU^ a \t) + dM^\t)\ (18) 
is the strain tensor at r^ a \ Since v\j > = we can rewrite the interaction Hamiltonian as 

iV 



flint = ~J2 [^fdiUj^^a^] . (19) 



a=l 



18 



The full Hamiltonian is now 



H 



d 3 r 



v 

N 



P.: 

2 



Uitii + -CijkidiUjd k ui 



N 



a=l 



2^0 a z + 2 b x 
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(20) 



Q = l 



For weak TLS-wave system interactions we can in principle start with this Hamiltonian and 
derive a master equation for the observed flexural wave mode of interest that interacts with 
the N TLS's. The rest of the elastic wave normal modes then form the TLS bath, as well as 
mediate the interactions between the TLS's. Instead, we will adopt a less rigorous approach 
to derive the approximate form of the elastic wave-induced interaction between any pair 
of TLS's. We assume that the timescale for the phonon mediated interaction between two 
TLS's is much shorter than their internal dynamics timescale; the two TLS's are therefore 
approximated as 'frozen,' with each in a given spin state. We take as our starting point the 
following Hamiltonian for the interaction between two TLS's without the non-interacting 
TLS part: 

2 



H; 



approx2TLS 



I d 3 r ^ + 

^ - " a=l 



(21) 



Next, we express this approximate Hamiltonian operator at t = in terms of the normal 
mode, phonon creation and annihilation operators. We define 



uAr 



and 




(22) 



(23) 



where [<Za,aL] = Sp,p>, with (3 labeling the normal mode. The normal modes are solutions 
to 

Cijkidjd k up t i = -pupupj. (24) 



Substituting Eqs. (|22|) and (|23|) into Eq. (21) and using Eq. (24) and the orthonormality 

dhup^ug./f) = S PiP > (25) 



and completeness relations 



v 
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and 

respectively, we obtain 

-f^approx2TLS 

where the function f« is defined as 



E 

/8 



hu 



13 (apa\ + a^ap) + fpap + j*p^p 



(26) 



(27) 



h 



(28) 



We now redefine the creation/annihilation operators so as to get rid of the linear operator 



terms in Eq. (27). We make the substitution dp = bp + cp, where we have included hats 



here to emphasize that a and b are operators, whereas c is a commuting number. We obtain 
(dropping the hats) 

hur- 



H. 



approx2TLS 



E 





+fpbp + ffll + cpfp + c* p f*p + hupcpcp 



(29) 



Defining cp = —f^/huj and c% = —fp/huj, we see that the linear operator terms in Eq. (29) 
cancel out and we have 



H, 



approx2TLS 



E 



(30) 



The TLS-TLS interaction term we are seeking is contained within the quadratic /-term in 



Eq. (|30|). Substituting in Eq. (|28|) and simplifying, we obtain for the TLS-TLS interaction: 

(31) 



#TLS-TLS = -O-^&VikVjl 



where we have neglected TLS self-interaction terms. 

The strength of the interaction between the two TLS's will depend on the nature of the 



elastic medium in which the TLS's are embedded, as expressed by the mode sum in Eq. (31 ). 



Let us now try to come up with a simple semiquantitative approximation to the mode sum 



part in the interaction term (31 ) using dimensional analysis. From the completeness relation 



(26), the displacement mode function upj has the dimensions L 3 / 2 in terms of some to-be- 



determined length scale L. The mode frequency depends on the speed of sound v, and so 
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scales as up ~ v/L. Thus, the overall length dimension for the mode sum in Eq. (31) is 
L~ 3 . The relevant length scale, however, depends on the geometry of the embedding elastic 
medium. For a bulk, three-dimensional (3D) medium where the two TLS's are far from any 
of the medium boundaries, the appropriate length scale must be the separation ryi between 
the two TLS's. Thus, for a 3D medium we have 

o-3D ^ (1) Jft]^_ J_ fW) 
^2TLS ~ z a z 2 r 3 ' 
P u '12 

where we have neglected the anisotropy of the deformation potential. For a membrane- 
like elastic medium, where the separation between the two TLS's is large compared to the 



membrane thickness d, we must lose one of the factors in (32), to be replaced by d. Thus, 
for an effectively 2D medium, we have 

7T2D ^(1)^(2) 1 / QQ \ 

H ™~ aM ^dr{ 2 (33) 

Finally, for a wire-like elastic medium where the separation between the two TLS's is large 
compared to the wire's crossectional dimensions d and w, we must lose two of the r 12 factors 



in (32). Thus, for an effectively ID medium, we have: 



^ms~^? , -S-— ■ (34) 

Note that, as the dimensions of the elastic structure are reduced, the TLS-TLS interaction 
becomes longer ranged. In particular, for a wire-like structure, the reduced volume and 
hence reduced number of TLS's will in part be compensated by a longer ranged interaction. 

B. Three Interacting TLS's 



We now include TLS-TLS interactions. We group all variables in Eq. (31) except the 
sigma operators into a single variable, ^ a ^\ simplifying the TLS-TLS interaction Hamilto- 
nian to 

^TLS-TLS = -otooP(W. (35) 

For the plots in this section and the next, we choose a value for ( and then generate a 
random within ±50% of this value for each pair of TLS's. Unless otherwise specified, 
the values are centered around ( = 0.1/6. We plot P n vs ut for a resonator coupled to three 
interacting TLS's. 
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Fig. 17 shows that P n decays exponentially as a function of time, and we again apply a 



linear fit to the log plot to extract a decay rate. In Fig. 18 we plot the normalized decay 
rate for a resonator coupled to three non-interacting TLS's (solid), to three interacting 
TLS's (dash), and to an Ohmic bath (dot-dash). The dot-dashed and dashed curves are 
practically indistinguishable, suggesting that the addition of TLS-TLS interactions allows 
the three TLS's to absorb energy like an Ohmic bath, even for higher-n Fock states. 
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FIG. 17: (Color online) Left: P n vs cut for a resonator coupled to three interacting TLS's. Right: 
Log of P n vs Lot. For all curves T\ = 10. 



We now investigate the decay of a superposition state. Fig. [19] shows the on- and off- 
diagonal decay times Tq h and 2T nn as a function of n for a range of TLS T\ values. As in 
the case of three non-interacting TLS's, the curves have slopes ~ —1 and also T 0n w 2T nn ; 
dephasing is negligible. 



Fig. 20 shows the T\ dependence of the on- and off-diagonal decay times for a resonator 



coupled to three non-interacting (black) and three interacting (gray) TLS's. The decay times 
are reduced for the resonator coupled to interacting TLS's, with the same unexpected T\ 
dependence as noted in the previous section. 

C. Six Interacting TLS's 



We now couple the resonator to six interacting TLS's. Fig. 21 shows the number state 
probability as a function of time for a resonator coupled to six damped, interacting TLS's. 
The shape of the curves is similar to that for six non-interacting TLS's, with the log plot 
appearing approximately linear. Again, the oscillations appearing in the larger n curves 
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11 



without interactions 




FIG. 18: (Color online) Normalized decay rate vs n for single Fock states. The resonator is 
coupled to three non-interacting TLS's (solid), three interacting TLS's (dash), and to an Ohmic 
bath (dot-dash). For all curves T\ = 10. 
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FIG. 19: (Color online) T$ n and T nn for a superposition state l^o) = l/v2(|0) + |n)) for a range of 
T\ values. The resonator is coupled to three interacting TLS's. For the final curve the TLS's are 
not damped, and thus a T\ time is not given. 
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FIG. 20: T nn and T$ n vs T\ for the n = 4 superposition state. The resonator is coupled to three 
non- interacting (black) and three interacting (gray) TLS's. 



at long times are numerical artifacts due to the exponentially small decay probabilities. In 



Fig. 22 we plot the decay rate as a function of n for six TLS's with (dash) and without (solid) 
TLS-TLS interactions, and for a resonator coupled only to an Ohmic bath (dot-dash). We 
note that the decay is similar for the two cases with slope close to one. 
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FIG. 21: (Color online) Left: P n vs uot for a resonator coupled to six interacting TLS's. Right: 
Log of P n vs u)t. For all curves T\ = 10. 

Next we study the T\ dependence of the on- and off-diagonal terms of the density matrix 



for a superposition of Fock states, as we did in Fig. 20. Fig. 23 shows T nn and T 0n as a 
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function of T\ for the n = 4 superposition state for a resonator coupled to six non-interacting 
(black) and six interacting (gray) TLS's. The plot shows that in both cases a reduction of 
T\ causes an increase in the decay time of the on- and off-diagonal terms. As for the three 
interacting TLS case, the addition of TLS-TLS interactions decreases the decay times. 
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FIG. 22: (Color online) Normalized decay rate vs n for a single Fock state. The resonator is 
coupled to six non-interacting TLS's (solid), six interacting TLS's (dash), and to an Ohmic bath 
(dot-dash). For all curves T\ = 10. 



Lastly, we plot T nn and T 0n as a function of the TLS-TLS coupling parameter (. Fig. [24 
shows the C dependence of the decay times for two different sets of random ( ( - a/3 \ In both 
cases the decay time of the diagonal terms, T nn , shows a linear dependence on the strength 
of the TLS-TLS coupling, with a slight variation in the slope for the two realizations. The 
off-diagonal terms, particularly for the second group of random (( a ^ values, decay less 
uniformly with respect to £, but the overall behavior shows a clear dependence on (, with 
stronger TLS-TLS coupling leading to faster decay of both the diagonal and off-diagonal 
terms of the density matrix. 



V. CONCLUSION 



In this work we have explored the effects of TLS's on damping and decoherence of Fock 
states in a mechanical resonator at low temperatures. We began our investigation with 
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FIG. 23: T nn and To„ vs T\ for the n = 4 superposition state. The resonator is coupled to six 
non-interacting (black) and six interacting (gray) TLS's. 

a resonator coupled to a single TLS, and then increased the number to three, and then 
six TLS's. For Fock states in a resonator coupled to a single, damped TLS, we observed 
amplitude-dependent damping, an indication of TLS saturation. For a resonator coupled to 
three TLS's we found that the decay of Fock states and Fock state superpositions was similar 
to that due to an Ohmic bath, particularly when TLS-TLS interactions were included. We 
noted that there was still some variation in the decay for different realizations of the random 
TLS variables, reflecting the fact that we were between a single TLS and a dense spectrum 
of TLS's. For a resonator in the presence of six TLS's we found that the damping of a Fock 
state went approximately as Ti/n, as expected for an Ohmic bath. Further, for the decay 
of a superposition of Fock states we found that the off-diagonal terms of the density matrix 
decayed twice as slowly as the diagonal terms, as expected for a resonator coupled to a bath 
of free oscillators. While the on- and off-diagonal terms showed an unexpected dependence 



on Ti, the analysis in Sec. |III D] suggested a possible explanation for their behavior in terms 
of T\ dependences of the TLS decay line widths. 

This work highlights the need for analytical approximations in order to understand the 
numerical results. One possibility is the use of a polaron-like transformation to take into ac- 
count the correlation between the TLS's and the resonator.^ This may allow us to simplify 
by approximation the equations of the oscillator and damped TLS's in order to understand, 
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FIG. 24: T nn and Tb n vs the TLS-TLS interaction strength £ for two different realizations of the 
random values. The resonator is initially in the n = 4 superposition state. For all curves 
Ti = 10. 

for example, the qualitative saturation dependence. Furthermore, the Ohmic bath-like de- 
pendence, i.e., decay rate proportional to n, for several near-resonant interacting TLS's 
suggests that appropriate analytical methods can help understand the relevant simulations. 
Much work remains to be done to understand the quantum-classical correspondence, par- 
ticularly now that experiments demonstrating mechanical systems in the quantum limit are 
becoming a reality. 
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